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Abstract. We use numerical simulations of galaxy clus- 
ters in different cosmologies to study their ability to form 
large arcs. The cosmological models are: Standard CDM 
(SCDM; = 1, f^A = 0); rCDM with reduced small-scale 
power (parameters as SCDM, but with a smaller shape 
parameter of the power spectrum); open CDM (OCDM; 

= 0.3, 17a = 0); and spatially flat, low-density CDM 
(ACDM; flo = 0.3, JIa = 0.7). AU models are normalised 
to the local number density of rich clusters. Simulating 
gravitational lensing by these clusters, we compute optical 
depths for the formation of large arcs. For large arcs with 
length-to-width ratio > 10, the optical depth is largest for 
OCDM. Relative to OCDM, the optical depth is lower by 
about an order of magnitude for ACDM, and by about two 
orders of magnitude for S/rCDM. These differences origi- 
nate from the different epochs of cluster formation across 
the cosmological models, and from the non-linearity of the 
strong lensing effect. We conclude that only the OCDM 
model can reproduce the observed arc abundance well, 
while the other models fail to do so by orders of magni- 
tude. 

Key words: Galaxies: clusters: general — Cosmology: 
dark matter — Cosmology: gravitational lensing — Cos- 
mology: large-scale structure of Universe 



1. Introduction 

Galaxy clusters form differently in different cosmological 
models. Their formation history and their internal struc- 
ture are influenced by the cosmological parameters. In 
dense model universes, fio ^ 1, clusters form at signifi- 
cantly lower redshifts than in low-density model universes 
(e.g. Richstone, Loeb, & Turner 1992; Bartelmann, Ehlers, 
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& Schneider 1993; Lacey & Cole 1993, 1994). The cosmo- 
logical constant has a fairly moderate influence on the 
formation timescale. Delayed formation is reflected in the 
abundance of cluster substructure. Moreover, central clus- 
ter densities are higher in clusters that form earlier. 

It is currently unclear whether the different degrees 
of cluster substructure expected in cosmological models 
with different mean densities lead to observational conse- 
quences that can signiflcantly distinguish between high- 
and low-density universes. While earlier studies found 
cluster X-ray morphologies and density profiles to differ 
significantly between different cosmologies (Evrard et al. 
1993; Crone, Evrard, & Richstone 1994; Mohr et al. 1995; 
Crone, Evrard, & Richstone 1996), more recent work con- 
cluded that X-ray morphologies of clusters at the present 
epoch are fairly similar in different cosmological mod- 
els, rendering significant distinctions difficult (Jing et al. 
1995). The issue of constraining cluster shapes using their 
X-ray emission was also addressed in detail by Buote & 
Tsai (1995a, b). The weak gravitational lensing effect al- 
lows a measurement of the morphology of the projected 
mass (Wilson, Cole, & Frenk 1996; Schneider & Bartel- 
mann 1997) and constrains cluster density profiles (Crone 
et al. 1997), but it remains to be shown with realistic nu- 
merical cluster models whether this method provides a 
more sensitive tool to quantify cluster morphologies than 
that provided by X-ray emission. 

An alternative tool is offered by the strong gravita- 
tional lens effect. In order to be strong lenses (i.e. in order 
to produce appreciable numbers of large arcs from back- 
ground sources), clusters have to satisfy several criteria. 
First, they need to be compact, that is, their central sur- 
face mass densities need to surmount the critical surface 
mass density for lensing. The latter depends on redshift. 
For background sources at redshifts Zg ^ 1, clusters at 
redshifts 0.2 ^ Zc ^ 0.4 are the most efficient lenses. If 
arcs are to be produced in abundance, a sufficiently large 
number of concentrated clusters must be in place at those 
redshifts. Second, strong lensing is a highly non-linear ef- 
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feet. This is mainly because the number of strong lens- 
ing events depends sensitively on the number of cusps 
in, and the length of, the caustic curves of the lenses. 
Cusps require asymmetric lenses. Asymmetric, substruc- 
tured clusters are thus much more efficient in producing 
large arcs than symmetric clusters, provided the individ- 
ual cluster sublumps are compact enough (Bartelmann, 
Steinmetz, & Weiss 1995). Both arguments show that the 
influence of cosmology on the structure and the formation 
timescale of clusters should strongly affect their ability to 
form large arcs. Clusters that are being assembled from 
compact subclusters at redshifts where lenses are most ef- 
ficient, 0.2 ^ z ^ 0.4, should produce many more arcs 
than clusters which form at later redshifts from sublumps 
which are less compact. 

Does this line of reasoning imply that different cosmo- 
logical models can be distinguished by the number of arcs 
that are expected in them? More precisely, does the higher 
compactness of clusters in low-density universes, and the 
later formation time of clusters in high-density universes, 
lead to such different numbers of large arcs that limits on 
cosmological parameters could be obtained from counting 
arcs? This is the question addressed by this paper. 

In order to pursue it, we use galaxy clusters simulated 
in a variety of cosmological models. The simulations are 
described in Sect. The simulated clusters are then in- 
vestigated as to their strong-lensing effects, following the 
prescription in Sect. ||. Results are presented in Sect. ^ 
and summarised in Sect. ^. The paper concludes with a 
discussion in Sect. ||. 

Wu & Mao (1996) already considered the influence of 
A on arc statistics, however with spherically symmetric, 
non-evolving clusters. They found that this cluster model 
predicts ~ 2 times more arcs in a low-density, spatially 
flat universe (fio = 0.3, f^A = 0.7) than in an Einstein-de 
Sitter universe, and that the latter model falls short by a 
factor of ~ 4 to explain the observed arc abundance. Using 
a similar model, Hamana & Futamase (1997) pointed out 
that the expected number of observed arcs increases when 
the evolution of the background sources is taken into ac- 
count. Hattori, Watanabe, & Yamashita (1997) included 
observational selection effects in a study of arc statistics 
using spherically symmetric cluster models. Using CDM 
cluster models at a single redshift in an Einstein-de Sit- 
ter universe, van Kampen (1996) studied the influence 
of the normalisation of the power spectrum on the num- 
ber of arcs per cluster and concluded that the normalisa- 
tion should be somewhat higher than derived from cluster 
abundance. 

2. Cluster simulations 

We use two different sets of cluster simulations. A first set 
of five clusters is taken from high-resolution cosmological 
simulations kindly made available by the GIF collabora- 
tion, and a second set of four clusters was simulated for 



the purpose of studying cluster mass profiles in different 
cosmologies. Different numerical techniques were used for 
the two simulation sets. They have in common that within 
each set, the same random phases are used for the initial 
density field in all cosmological models studied, so that 
the clusters can be compared individually. The normalisa- 
tion of the power spectra agrees approximately with the 
normalisation to the local cluster number density. Apart 
from improved statistics, the approach using two differ- 
ently simulated cluster sets allows us to test whether dif- 
ferent numerical techniques yield different results for arc 
statistics. Although the normalisations for the two simu- 
lation sets are slightly different, the arc statistics obtained 
from the two sets individually agree well with each other. 
The clusters show density profiles pir) that are well fit- 
ted by the two-parameter function suggested by Navarro, 
Frenk, & White (1996), 



Navarro et al. (1996) introduced the concentration param- 
eter c = 7*200 /''s, where r2oo is the radius enclosing an 
average overdensity of 200 times the cosmic background 
density. At 2 = 0, we find 

r 5 for S/rCDM 
c w <^ 7 for ACDM . (2) 
[9 for OCDM 

The different values for c reflect the different cluster for- 
mation times. When clusters form earlier, their concentra- 
tion is higher. 

2.1. First cluster sample: GIF simulations 

2.1.1. The GIF project 

The GIF project is a joint effort of astrophysicists from 
Germany and Israel. Its primary goal is to study the for- 
mation and evolution of galaxies in a cosmological context 
using semi-analytical galaxy formation models embedded 
in large high-resolution iV-body simulations. This is done 
by constructing merger trees of particle haloes from dark- 
matter only simulations and placing galaxies into them 
using a phenomenological modelling (for a detailed de- 
scription of this procedure as well as results cf. Kauff- 
mann et al. 1997). In order to achieve both good statistics 
and an accurate treatment of early epochs, high resolu- 
tion simulations are needed which nevertheless contain a 
fair sample of the Universe, thus accounting correctly for 
the influence of large-scale structure on galaxy formation. 
Those characteristics also make these simulations suitable 
for the present project. 

2.1.2. The simulations 

The code used for the GIF simulations is called Hydra. It is 
a parallel adaptive particle-particle particle- mesh (AP^M) 
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code (for details on the code cf. Couchman, Thomas, & 
Pearce 1995; Pearce & Couchman 1997). The current ver- 
sion was developed as part of the VIRGO supercomputing 
project and was kindly made available by them for the 
GIF project. The simulations were started on the CRAY 
T3D at the Computer Centre of the Max-Planck Society 
in Garching (RZG) on 128 processors. Once the clustering 
strength required an even larger amount of total memory, 
they were transferred to the T3D at the Edinburgh Par- 
allel Computer Centre (EPCC) and finished on 256 pro- 
cessors. 

A set of four simulations with N = 256'^ and with 
different cosmological parameters was run. Apart from 
the fiducial Cold Dark Matter (CDM)^cenario, denoted 
SCDM, which has ^ 1 and h = 0.50, another f^o = 1 
and h = 0.5 model was run (rCDM) which has the same 
shape parameter for the power spectrum, P = 0.21, as 
the remaining two models. Those are both models with 
r^o = 0.3, the first one being a flat model with a cosmo- 
logical constant (ACDM, JIa = 0.7, h = 0.7), and the last 
model being an open model (OCDM, = 0, h = 0.7). 
The case of the rCDM model is particularly interesting 
because it shares the ^Iq = 1 dynamics with the SCDM 
model, but has a power spectrum with the same shape 
parameter as the two low-ilp models. A value of P = 0.21 
is usually preferred by analyses of galaxy clustering, cf. 
Peacock & Dodds (1994). This is achieved in the rCDM 
model despite ilo = 1 and h = 0.5 by assuming that a 
massive neutrino (usually taken to be the r neutrino) had 
existed during the very early evolution of the Universe. 
It must have decayed later, thus shifting the epoch when 
matter started to dominate over radiation in the Universe, 
and the neutrino mass and lifetime are chosen such that 
P = 0.21. For a detailed description of such a model see 
White, Gelmini, & Silk (1995). 

The GIF simulations adopt the power spectrum 
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6.4P^i/i"iMpc 
3.0P"^/i"^Mpc 

1.13 
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(Bond & Efstathiou 1984). In order to completely fix 
P{k), the normalisation. A, has to be chosen. This can 
be done on the basis of measurements of the microwave 
background anisotropics by the COBE satellite. However, 
this approach suffers from the fact that COBE measured 



fluctuations on scales much larger than those pertinent to 
the simulations. Hence, one has to assume that there is 
no additional physics that could alter the result like, e.g., 
gravitational waves or a slight tilt of the initial spectrum 
away from the scale-invariant form. The approach taken 
here avoids this problem. The mass function of objects in 
the Universe is very steep at the high-mass end. In other 
words, massive objects (like clusters of galaxies) are not 
only rare, but their abundance sensitively depends on the 
amplitude of the power spectrum. White, Efstathiou, & 
Frenk (1993) introduced this way of fixing the amplitude 
by determining erg, the square root of the variance of the 
density field smoothed over 8/i~^Mpc spheres, such that 
the observed abundance of rich clusters is matched. They 
used the cluster mass function. Recent studies of the clus- 
ter X-ray temperature function (Eke, Cole, & Frenk 1996; 
Viana & Liddle 1996) find similar results. For the low- 
density GIF simulations, the result by Eke et al. (1996) 
was taken, 



(0.52 ±0.04) 17^0.46+0. 10 Oo 

(0.52 ±0.04) f]-0-52+0.13 0o 



for r^A 

for ilo - 



= 

f^A 



(5) 



For the ^Iq — 1 simulations, slightly larger values than 
suggested by eq. (||) were adopted, according to the ear- 
lier result by White et al. (1993). Table summarises the 
model parameters. 



Table 1. Cosmological parameters of the GIF models, ilo and 
S7a are the density parameters for matter and the cosmologi- 
cal constant, h is the Hubble parameter, as is the variance of 
the density field in spheres of 8 Mpc, and F is the shape 
parameter of the power spectrum. Also given are the size of 
the cosmological simulation box and the mass Mmax within 
1.5 Mpc radius of the most massive cluster in 10^^ Mq. 



Model 




Ha 


h 


as 


r 


Box Size 
















[Mpc/h] 




SCDMl 


1.0 


0.0 


0.5 


0.60 


0.50 


85 


0.74 


rCDMl 


1.0 


0.0 


0.5 


0.60 


0.21 


85 


0.74 


ACDMl 


0.3 


0.7 


0.7 


0.90 


0.21 


141 


0.84 


OCDMl 


0.3 


0.0 


0.7 


0.85 


0.21 


141 


0.85 



As usual, the Hubble constant is written as Ho — 100 h 
kms"^ Mpc"\ 



The parameters shown in Table |i| were chosen not 
only to fulfil cosmological constraints, but also to al- 
low a detailed study of the clustering properties at very 
early redshifts. The masses of individual particles are 
1.0 X lOio/i-iAfo and 1.4 x 10^" Mq for the high- 
and low-f^o models, respectively. The gravitational soft- 
ening was taken to be 30 kpc. 

Clusters are obtained from the simulation as fol- 
lows. High-density regions are searched using a standard 
friends-of-friends group finder with a linking length of 
b = 0.05 times the mean interparticle separation. This se- 
lects only the dense cores of any collapsed object. Around 
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the centres of these, all particles are collected which lie 
within a sphere of radius ta = 1.5 ft." ^ Mpc, which corre- 
sponds to Abell's radius. These objects are taken as clus- 
ters. For our analysis, the five most massive clusters are 
cut out of the simulation volumes. This procedure needs to 
be expanded if the centres of two large clusters are closer 
together than r^^. In this case usually the more massive 
cluster is taken, and the other one is deleted from the list. 
In our case, however, this problem did not occur. The ini- 
tial density fields for the different cosmologies share the 
same random phases. 

2.2. Second cluster sample 

The second cluster sample is simulated using a special 
multi-mass technique which is explained in detail in Huss, 
Jain, & Steinmetz (1997). In contrast to the GIF simu- 
lations, this technique gives only one massive cluster per 
run. However, it allows one to study the evolution of one 
individual cluster without the need for extensive computer 
resources. 

The essential part of the multi-mass technique is the 
initial particle arrangement. It consists of three spherical 
layers embedded into a cubic volume, each filled with par- 
ticles of different mass. The central sphere encompasses 
the least massive particles and is surrounded by two shells 
of more massive particles. The rest of the cubic simula- 
tion volume is filled up with the most massive particles. 
The inner sphere must initially be large enough to enclose 
all particles which end up in a cluster. The gravitational 
forces on the particles are calculated using a combined 
GRAPE/PM A^-body code assuming periodic boundary 
conditions. The PM part performs force calculations with 
periodic boundary conditions for all particles. In the in- 
ner three shells, the force is additionally calculated with 
a PM code using vacuum boundary conditions. This force 
is subtracted from the periodically computed PM force to 
obtain the periodic contribution to the force only. This is 
added to the highly resolved force provided by the GRAPE 
board for the particles in the inner shells. 

The second cluster sample consists of 12 clusters in 
total. Four clusters are simulated for each of three differ- 
ent cosmologies, which resemble the SCDM, ACDM, and 
OCDM models of the GIF project. The model parame- 
ter are summarised in Table ||. All clusters belonging to 
one cosmological model are part of the same realisation of 
the corresponding density field. In addition, the phases of 
the initial Gaussian random field are identical in the three 
cosmological models. 

The initial conditions are calculated using eq. (|^) to 
match the power spectrum of the different cosmologies. 
The normalisation of P{k) is chosen as determined by 
White et al. (1993) for SCDM2, and for 0CDM2 and 
ACDM2 it matches the COBE CMB anisotropy measure- 
ments. Contributions from gravitational waves need not be 
considered since eq. (p|) assumes that the primordial power 



Table 2. Cosmological parameters of the models for the second 
cluster sample. The meaning of the symbols is the same as in 
Tab. |l[ 



Model 


no 


Ha 


h 


OS 


r 


Box size 
















[Mpc//i] 




SCDM2 


1.0 


0.0 


0.5 


0.60 


0.50 


144 


0.75 


ACDM2 


0.3 


0.7 


0.7 


1.12 


0.21 


201.6 


1.67 


OCDM2 


0.3 


0.0 


0.7 


1.12 


0.21 


201.6 


1.41 



spectrum after inflation is of the Harrison-Zel'dovich form. 
The COBE normalisation yields a slightly higher (Tg than 
that derived by Eke et al. (1996). The box size L of the 
simulation volume is fixed to 288 Mpc in physical coor- 
dinates for each run to achieve the same physical spatial 
resolution in all models. The mass resolution in the central 
sphere is 4.9 x 10^" Mq for the high-density model, 
and 2.9 x 10^" h^^ Mq for the low-density models. 

Since the simulations result in one massive cluster in 
the high-resolution region, it can easily be identified by 
looking for the deepest potential well in this region. How- 
ever, particular attention has to be paid in order to set 
up initial conditions with a suitable overdense region in 
the central sphere, representing the seed for the massive 
cluster. This is done in the following way. First, a pre- 
simulation is performed by filling the whole simulation 
box with particles of the same mass as those in the sec- 
ond layer. At z = 0, the cluster-like objects are identified 
using a special group-finding algorithm (Huss et al. 1997). 
The particles finally forming these objects define the corre- 
sponding overdense region at the initial redshift. The final 
starting configuration is then centred on one of these re- 
gions. By adding small scale power to such a density peak, 
the clustering properties of the region can be changed. 
Hence, it is possible that the simulation finally arrives at 
several low-mass objects rather than at one massive halo. 
This can be avoided by testing the clustering properties 
using the Zel'dovich approximation. When propagated to 
low redshifts with this technique, the particles in the cen- 
tral sphere must form a distinct matter accumulation in 
the centre rather than showing only filamentary structure. 

With this procedure, four suitable overdense regions 
are identified in the SCDM2 model. For the ACDM2 and 
for the 0CDM2 clusters, the starting configurations are 
centred on the same regions as for the SCDM2 clusters. 
This is possible since the clustering properties are defined 
mainly by the local realisation of the random field, which 
is the same for all three models. However, the final cluster 
haloes need not represent the most massive clusters in the 
simulation box. 

3. Simulations of arcs 

Our method to investigate the arc-formation statistics of 
the numerical cluster models was described in detail by 
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Bartelmann & Weiss (1994). We will therefore keep the 
present description brief and refer the reader to that paper 
for further information. For general information on gravi- 
tational lensing, see Schneider, Ehlers, & Falco (1992) or 
Narayan & Bartelmann (1997), and references therein. 

The numerical cluster models yield the spatial coordi- 
nates and velocities of discrete particles with equal mass. 
In order to use them for gravitational lensing, we need 
to compute the surface mass density distribution of each 
cluster model in each of the three independent directions 
of projection. The mass density is first determined by sort- 
ing the particles into a three-dimensional grid and subse- 
quently smoothing with a Gaussian filter function. The 
grid resolution and the width of the Gaussian are adapted 
to the numerical resolution of the A'^-body codes in or- 
der not to lose spatial resolution by the smoothing of the 
density field. The smoothed density field is then projected 
onto the three sides of the computation volume to obtain 
three surface-density fields for each cluster. 

The physical surface mass density fields are then scaled 
by the critical surface mass density for lensing, which 
apart from the cosmological parameters depends on the 
cluster- and source redshifts. We keep the redshift for all 
sources fixed at Zg = 1, and the clusters are at < Zc < 
1. This finally yields three two-dimensional convergence 
fields k{x, Zc) for each cluster model at redshift Zc- From 
k(cc, Zc), all quantities determining the local lens mapping, 
i.e., the deflection angle and its spatial derivatives, can be 
computed. We determine the lens properties of the clus- 
ters on grids with an angular resolution of 0'.'3 in the lens 
plane in order to ensure that lensed images be properly 
resolved. 

Sources are then distributed on a regular grid in the 
source plane. The resolution of this source grid can be kept 
low close to the field boundaries because there no large 
arcs occur. Close to the caustics of the clusters, where 
the large arcs are formed, the source-grid resolution is in- 
creased with the increasing strength of the lens. For our 
later purpose of statistics, sources are weighted with the 
inverse resolution of the grid on which they are placed. 
The sources are taken to be intrinsically randomly ori- 
ented ellipses with their axis ratios drawn randomly from 
the interval [0.5, 1], and their axes determined such that 
their area equals that of circles with radius 0'.'5. Although 
this choice of source properties appears fairly simple, it 
should not affect the arc statistics because these mainly 
reflect the local properties of the lens mapping, which are 
independent of the particular choice of source size or ellip- 
ticity distribution. We checked that a change in average 
source size did not change the results. 

The sources are then viewed through the cluster lenses. 
All images are classified in the way detailed by Bartel- 
mann & Weiss (1994). Among other things, the classifica- 
tion yields for each image its length L, its width W, and 
its curvature radius R. In total, we classify the images of 
about 1.3 X 10^ sources. 



Knowing the area covered by the cluster fields, and 
having determined the frequency of occurrence of image 
properties such as a given length, width, and curvature 
radius, we can compute cross sections a for the formation 
of images with such properties. The arc cross section of a 
cluster is defined as the area in the source plane within 
which a source has to lie in order to be imaged as an arc. 
We mostly focus on cross sections for the length-to-width 
ratio r of arcs. Apart from the image properties, the cross 
sections depend on redshift, a = a{z). 

Given cross sections cr(z), we compute the optical 
depth r for the formation of large arcs. The optical depth 
is the fraction of the entire source plane which is covered 
by cluster cross sections, 

nc r' , ^^ AV{z) , , 

where ric is the present cluster number density, Ds is the 
angular-diameter distance to the source plane, and dV^(z) 
is the proper volume of a spherical shell of width dz about 
z. The factor (1 + zY accounts for the cosmological ex- 
pansion factor. 

4. Results 

J^.l. Cross sections 

According to the prescription in Sect. ||, we first calculate 
cross sections for each individual model cluster, projected 
in each of the three independent spatial directions. Inter- 
polating in redshift between the redshifts of the model 
clusters, this yields cross sections cr(z) as a function of 
redshift. We then average these cross sections (1) over the 
three projection directions and (2) over all model clus- 
ters within a given cosmological model. We thus obtain 
cross sections for the four cosmological models. Figure |^ 
shows an example, the cross sections for arcs with length- 
to-width ratio r > 7.5. 

The averaged cross sections in Fig. |l] reveal huge dif- 
ferences between the cosmological models. While the cross 
sections for standard CDM (SCDM) and rCDM are com- 
parable, the maximum cross sections for ACDM and open 
CDM (OCDM) exceed that for SCDM by about half 
and one order of magnitude, respectively, and the red- 
shift range where (t(z) > is wider in 0/ACDM than 
in S/rCDM. Cross sections for other, large values of the 
length-to-width ratio r, or for large arc lengths, show a 
qualitatively similar behaviour. 

Since the clusters in different cosmologies arise from 
initial density perturbations with the same random 
phases, they can also be compared individually rather 
than statistically. On the whole, the individual clusters 
show the same qualitative behaviour as the averaged cross 
sections shown in Fig. |l]. Results obtained for each cluster 
set individually agree well with each other. 
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Fig. 1. Averaged cross sections for arc length-to- width ratio 
r > 7.5 for clusters in the four different cosmological models, 
distinguished by line types as indicated. The figure shows that 
clusters in the rCDM model produce the fewest arcs, and clus- 
ters in the open CDM model the most. Note the logarithmic 
scale of the ordinate: The maxima of the cross sections differ by 
more than an order of magnitude. The redshift ranges where 
a(z) > are larger in O/ACDM than in S/tCDM. 



4-. 2. Arc-cluster redshift 

At what redshifts do we expect to find the most clus- 
ters that produce large arcs? In other words, clusters at 
which redshift contribute most to the arc optical depth? 
To answer this question, we compute the optical depth 
T from eq. (^) and the differential optical depth dr/dz, 
and plot in Fig. ^the normalised differential optical depth 
dr/dz as a function of redshift for the four cosmolog- 
ical models. 

The curves in Fig. ^ show that the differential optical 
depth peaks around z ^ 0.3 — 0.4. The bars inserted in the 
figure indicate the 1-ct redshift range, and the dots show 
the average arc-cluster redshift. Although there is a slight 
tendency that the mean arc-cluster redshift is smallest in 
the SCDM model, larger for ACDM, and largest for rCDM 
and OCDM, the redshift variances are large enough for the 
redshift ranges in the cosmological models to overlap. The 
figure furthermore suggests that the differences between 
cosmological models are dominated by noise. Numbers are 
given in Table ||. 



Table 3. Mean redshifts and redshift ranges for clusters pro- 
ducing large arcs in the four different cosmological models. 



model 


Zc 




SCDM 


0.29 


0.09 


rCDM 


0.38 


0.06 


ACDM 


0.36 


0.09 


OCDM 


0.39 


0.12 



Fig. 2. Normalised, differential optical depth as a function of 
redshift, for arcs with length-to- width ratio r > 7.5. The curves 
indicate the most probable redshift for a cluster forming arcs 
for the four cosmological models used. The bars show the 1-a 
redshift range, the dots indicate the mean arc-cluster redshift. 
The plot shows that there is no significant difference in arc 
cluster redshift between the four cosmological models. 



4-3. Optical depth 

We now compare the optical depth r for formation of arcs 
with given length-to-width ratio r in the four cosmologi- 
cal models. As before, the optical depth is calculated from 
eq. (^). We do not specify the cluster number density ric 
yet, but calculate the optical depth per unit cluster den- 
sity, n~^T. Results are shown in Fig. || 

Figure ^ confirms the trends indicated by the cross 
sections in Fig. ^ but allows to compare optical depths 
for a wide range of arc length-to-width ratios r. There 
is an interval at intermediate r, 5 r < 10, where the 
optical depths for SCDM and tCDM are almost equal. 
Only at r <; 10 does the optical depth in tCDM models 
drop below that of SCDM models. For r ^ 4, the optical 
depth for ACDM models is constantly higher than that for 
SCDM models by a factor of ~ 10. The optical depth for 
the OCDM model is highest, exceeding the SCDM value 
by up to ^ 2 orders of magnitude at large r. 

The hatched region around the SCDM curve illustrates 
1-a bootstrap errors, which we obtained by bootstrapping 
the cluster sample. They give an impression of the uncer- 
tainty of the optical depth due to the limited number of 
clusters in our samples. The uncertainty due to the spe- 
cific realisation of the lensed background galaxy sample 
are smaller by about a factor of five. 

In order to emphasise the results. Fig. |^ shows the op- 
tical depths for tCDM, ACDM, and OCDM, divided by 
the optical depth for SCDM. 

5. Discussion 
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Fig. 3. Optical depth for arc formation as a function of arc 
length-to- width ratio r, for the four different cosmological mod- 
els. The optical depth for the rCDM and the standard CDM 
models are comparable for intermediate length-to-width ratios 
r. For large r, the optical depth is smallest for the rCDM 
model. The largest optical depth is produced by clusters in the 
open CDM model, followed by those in the ACDM model. For 
r ~ 10, the optical depths for open CDM, ACDM, and stan- 
dard CDM differ by about an order of magnitude each. The 
hatched area around the SCDM curve indicates 1-cr bootstrap 
errors. 

5.1. Results 

We have used numerical simulations to calculate the op- 
tical depth for the formation of large arcs in different cos- 
mologies. The simulated clusters are cut out of large cos- 
mological simulation volumes, so that the tidal effects of 
surrounding matter are taken into account. The cosmo- 
logical simulations were normalised to the observed lo- 
cal number density of rich galaxy clusters. Four different 
cosmogonic models were used. These are: standard CDM 
(SCDM), with no = 1, Qa = 0, h = 0.5, = 0.6, and 
shape parameter F = 0.5; a CDM model with reduced 
small-scale power (rCDM), which differs from SCDM only 
by the shape parameter F = 0.21: open CDM (OCDM) 
with n = 0.3, = 0, h = 0.7, cts = 0.9 or 1.1, and 
F — 0.21; and finally a spatially flat, low-density CDM 
model (ACDM) with Oq = 0.3, = 0.7, h = 0.7, 
as = 0.9 or 1.1, and F = 0.21. All cosmological simu- 
lations start from density perturbations with the same 
random phases, so that all clusters can be compared in- 
dividually in different cosmogonies. For SCDM, OCDM, 
and ACDM, we simulated nine clusters, and five clusters 
for rCDM. 

The lensing properties of the clusters with respect to 
large arcs were calculated in a way that has extensively 
been described earlier (Bartelmann & Weiss 1994; Bartel- 
mann et al. 1995). The calculations result in averaged 
cross sections for the clusters as a function of redshift, 
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Fig. 4. Optical depths for largo arcs, normalised by the arc 
optical depth in the standard CDM model. The curves empha- 
sise the large diflference in optical depth for the four models. 
For the largest r plotted, there arc more than two orders of 
magnitude between open CDM and standard CDM clusters, 
about one order of magnitude between ACDM and standard 
CDM, and clusters in the rCDM model are less efficient than 
the standard CDM models for r ^ 10. 



which can then be converted to optical depths t for the 
formation of large arcs. 

Our main result is that the optical depths for large 
arcs, with length-to-width ratio r > 10, differ by orders of 
magnitude for the different cosmologies. Generally, clus- 
ters in the SCDM and rCDM models produce the smallest 
optical depth. For r ^ 10, the optical depths for these two 
models are comparable, but for larger r, the optical depth 
in the rCDM model falls below that in the SCDM model. 
In the ACDM model, the arc optical depth is larger by 
about an order of magnitude than for SCDM, and the 
optical depth is largest in the OCDM model, exceeding 
the SCDM optical depth by about two orders of magni- 
tude. We emphasise that these results are independent of 
whether our cluster samples are in any sense complete or 
not, because the simulations are designed such that the 
clusters can be compared individually. Our conclusions 
do, however, rest on the assumption that the simulated 
clusters are typical for the clusters with the largest mass 
in each of the cosmological models. We believe that this is 
guaranteed by the large size of the cosmological simulation 
volumes from which the clusters were taken. 

It is a combination of effects that leads to the large dif- 
ference in arc optical depth across the cosmological mod- 
els that we have investigated, (i) Clusters form earlier in 
low-density than in high-density universes. In SCDM, nor- 
malised to the cluster abundance, the formation of such 
clusters which would in principle be massive enough for 
strong lensing is delayed to such low redshifts that they 
fail to be efficient lenses for sources at redshifts ~ 1- 
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(ii) For low-density universes, the proper volume per unit 
redshift is larger than for high-density universes. Given 
the observed number density of clusters today, this vol- 
ume effect increases the number of clusters between the 
source sphere and the observer when flo is small, (iii) Clus- 
ters that form early are more concentrated than clusters 
that form late. Clusters in low-density universes therefore 
reach higher central surface mass densities than clusters 
in high-density universes, (iv) Strong gravitational lensing 
is a highly non-linear effect. This is because the arc cross 
section of a cluster sensitively depends on the length of the 
caustic curve and the number of cusp points contained in 
it. The properties of the caustic curve do not only depend 
on the surface mass density, but also on the tidal field of 
a cluster, which is influenced by the cluster morphology, 
(v) Because of (iv), asymmetric clusters are much more 
efficient in producing large arcs than symmetric clusters 
(Bartelmann et al. 1995). The degree of substructure of a 
cluster is therefore very important for arc statistics. While 
clusters are in the process of formation, they are expected 
to be highly asymmetric. If this happens at redshifts where 
lensing is efficient for a given source population, the asym- 
metric cluster morphology further increases the strong- 
lensing cross section. Most clusters in OCDM and ACDM 
form at z ~ 0.3, exactly where lensing is most efficient 
for sources at ~ 1- Clusters in SCDM and rCDM form 
later, at z ~ 0.1, where their lensing efficiency and that of 
their sublumps is already suppressed by the lensing ge- 
ometry. Although clusters in S/rCDM form later than 
in 0/ACDM and should therefore be more asymmetric, 
the sublumps in O /ACDM clusters are more compact and 
thus tend to persist for a longer time after merging with 
the cluster. 



because the cross section per cluster should scale approx- 
imately with D^ff{z, Zs). This quantity is plotted in Fig. || 
for rio = 1 and = 0.3, with J^a = 0. Of course, this 
simple estimate completely neglects the influence of the 
change in cluster concentration across the cosmological 
models, and the non-linearities of the strong lensing effect. 
However, it suffices to demonstrate that large differences 
in the arc cross section are expected between high- and 
low-density universes. 




Fig. 5. Estimate for the number of efficient lensing clusters 
per redshift interval, dA^'ions/dz, as given in eq. (^). Results 
for fio = 1 and flo — 0.3, both with Ha = 0, are plotted, as 
indicated. The figure illustrates that the delayed formation of 
clusters in a high-density universe, combined with the effects 
of lensing efficiency and cosmic volume, already account for a 
large difference in the expected number of arcs. 



5.2. Illustration 



A simple Press-Schechter type argument illustrates the 
influence of formation time and cosmic volume. According 
to Press & Schechter (1974), the (comoving) fraction of the 
cosmic matter that is contained in clusters is 



F,{z) = - eric ( 

2 \V2anV{z) 



(7) 



with 6c ~ 1.686, (tr the variance of the density contrast on 
cluster scales today, and 'D{z) the (cosmology-dependent) 
linear growth factor of density perturbations. The cluster 
fraction at redshift z, normalised to the present cluster 
fraction, provides an estimate for the change in cluster 
number density with redshift. Multiplying with the proper 
cosmic volume An D^{z) |d(ci)/dz| dz of a shell of width dz 
and the squared effective lensing distance Dcs{z, Zg) yields 
an estimate for the number of efficient lensing clusters per 
redshift interval, 



diV, 



lens 



dz 



Fciz)xil+zfxD^^ff{z,z,)x47TD\z) 



diet) 



dz 



5.3. Influence of "missing clusters" 

All cosmological simulation volumes from which we have 
taken the cluster models are of order a few times 
10^ Mpc"^. Since the cluster mass function is very 
steep, we are therefore likely to miss the most massive 
clusters. In order to estimate their influence on the arc 
cross sections, we have repeated the arc simulations for 
SCDM with surface-mass densities rescaled to higher clus- 
ter mass. Let Mq and M > Mo be the original and rescaled 
cluster masses, respectively. Then, we take 



k' iX, Zn, M) ~ ( — 

^ ' ' ' * Mo 



,(8) 



1/3 



MoV^' 

X, Zr, Mn 

M j ' " 



(9) 



for the rescaled surface-mass density and compute the arc 
cross section from that. In effect, we calculate the arc cross 
section uizcTM) of a cluster of similar structure as the 
original one, but with higher total mass. For each cluster 
model at redshift Zc, we then average the arc cross sections 
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over mass, weighted by the Press-Schechter cluster mass 
f miction n(M), 



JdMn{M)c7{z^,M) 
/ dMn{M) 



(10) 



Finally, we compute mass-averaged optical depths (r) by 
substituting (cr(zc)) for (j{zc) in eq. (^). 

We found that (r) differs from r only by 10-15 per cent. 
Although the arc cross section is a fairly steep function of 
mass, the increase in a is more than compensated by the 
steep decrease of the cluster mass function. For example, 
quadrupling the masses of the rCDM clusters increases 
the averaged cross section for arcs with r > 10 by about 
two orders of magnitude, but decreases the cluster mass 
function by about three orders of magnitude, almost com- 
pletely cancelling the effect of the larger arc cross section. 
For clusters in the other cosmological models, the change 
in the optical depth should be smaller than that in rCDM. 
Clusters with higher masses than those contained in our 
sample can therefore safely be neglected. 

We have only studied the most massive clusters found 
in each cosmological simulation. Even then, the least mas- 
sive clusters in each sample contribute little or nothing to 
the arc optical depth. Extending our samples by includ- 
ing less massive clusters would therefore change the arc 
optical depth negligibly or not at all. 

Since we select the most massive clusters at redshift 
zero and study their progenitors, it could be that at 
higher redshift other clusters in the cosmological simu- 
lations would be more massive. In other words, it is not 
immediately clear that the progenitors of the most mas- 
sive clusters are also the most massive clusters present 
at higher redshift. In order to test that, we selected the 
rCDM model, where clusters form at the lowest redshifts, 
and checked whether the progenitors of the our sample of 
the five most massive clusters is identical with the sample 
of the five most massive clusters at the redshifts relevant 
for lensing. This turned out to be the case. In the other 
models, where clusters form at higher redshifts, possible 
misidentifications are even less likely than in rCDM. 

Given these results, we are confident that our cluster 
samples fairly reflect those clusters that dominate the op- 
tical depth for the formation of large arcs. 

6. Comparison with observations 

So can we constrain cosmological parameters through arc 
statistics? For that, we would have to compare the num- 
ber of observed arcs to that predicted by our models. 
This comparison is hampered by the fact that there is 
no complete sample of observed clusters selected by mass , 
as it should be for a fair comparison. There is one clus- 
ter sample, however, whose definition comes close to this 
criterion, namely the EMSS sample of X-ray bright clus- 
ters, for which the X-ray luminosity in the EMSS en- 
ergy band is Lx > 2 x 10'''* erg s"^ (h = 0.5, qo = 0.5). 



The number density of such clusters is estimated to be 
He ~ 2 X 10"^ /i^ Mpc"^ (Le Fevre et al. 1994). Arc sur- 
veys in this sample have shown that the number of arcs 
with r > 10 and a limiting magnitude of i? = 22.5 (or 
R — 21.5; these are the arc criteria set up by Wu & Ham- 
mer 1993) is roughly ~ 0.2 — 0.3 per cluster (Le Fevre et 
al. 1994; Gioia & Luppino 1994). 

Clusters with Lx > 2 x 10'*'' ergs" should be fairly 
represented by the massive simulated clusters in our sam- 
ples. Having velocity dispersions ^ 800kms~*, the em- 
pirical relation between velocity dispersion and X-ray lu- 
minosity obtained by Quintana & Melnick (1982) implies 
X-ray luminosities in the right range. We can therefore 
assume that the arc cross sections of our simulated clus- 
ters are typical for X-ray luminous clusters in the EMSS 
survey. 

The curves in Fig. ^give r. Using the number den- 
sity of bright EMSS clusters given above, 



-(r > 10) 



'2.9 X 10-^ (OCDM) 

3.3 X 10-^ (ACDM) 

4.4 X 10-^ (S/tCDM) 



(11) 



Since the whole sky has ^ 4.1 x 10^ square degrees, the 
total solid angle in which sources at Zg ~ 1 are imaged as 
large arcs with r > 10 is 



5lu 



1.2 X 10-1 sq. deg. (OCDM) 
1.4 X 10-2 sq. deg. (ACDM) 
1.8 X 10-3 sq. deg. (S/rCDM) 



(12) 



The sources which are imaged as arcs with the above prop- 
erties, r > 10 and R < 21.5 correspond to sources with 
R ^ 23.5 because of the magnification. Taking the num- 
ber densities compiled and measured by Small et al. , there 
are ^ 2 x lO'* such sources per square degree. The aver- 
age redshift of such sources is ~ 0.8 — 1 (e.g. Lilly et al. 
1995). Since the average arc-cluster redshift in our mod- 
els is at Zc ~ 0.3 — 0.4, the exact redshift of sources at 
z ~ 0.8—1 has only very little influence; the critical surface 
mass density changes by 10% when sources are shifted 
from z ~ 0.8 to z = 1.2. It follows that the number of 
such arcs on the whole sky expected from our simulations 
is 



r2400 (OCDM) 
iVarcs - < 280 (ACDM) 
[36 (S/rCDM) 



(13) 



There are ~ 7500 clusters on the sky which match the cri- 
teria of the EMSS bright cluster sample (Le Fevre et al. 
1994). Taking the number of arcs per cluster found in the 
EMSS clusters, the expected number of arcs on the whole 
sky is ^ 1500 — 2300. Despite the obvious uncertainties 
in this estimate, the only of our cosmological models for 
which the expected number of arcs comes near the observed 
number is the open CDM model. The others fail by one or 
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two orders of magnitude. The large differences in arc op- 
tical depth between the cosmological models investigated 
makes this result fairly insensitive to moderate uncertain- 
ties. It therefore appears fair to conclude that arc statistics 
in the framework of cluster-normalised CDM models de- 
mands that Qq is low, and that is small. Conversely, 
if r^o ^ 1, clusters have to form earlier than in our mod- 
els. We estimate with the simple Press-Schechter approach 
sketched above that in order to achieve that, CTg ~ 1.2 — 1.3 
would be necessary in the SCDM case. 

We have neglected the potential influence of cD galax- 
ies or cooling flows on the arc cross sections that could 
increase the central surface mass densities of the clusters 
and thus also their arc cross sections. Most probably, this 
influence is small compared to the huge differences be- 
tween the cosmological models. Nonetheless, we will study 
this issue in detail in a further paper because individual 
cluster galaxies may well affect arc morphologies, if not 
their total number. 
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